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Momentum Space 

Application to the Lippmann-Schwinger Resonating- Group Method and 

the pd Elastic Scattering 

Yoshikazu Fujiwara and Kenji Fukukawa 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

A practical method to solve cut-off Coulomb problems of two-cluster systems in the 
momentum space is given. When a sharply cut-off Coulomb force with a cut-off radius p is 
introduced at the level of constituent particles, two-cluster direct potential of the Coulomb 
force becomes in general a local screened Coulomb potential. The asymptotic Hamiltonian 
yields two types of asymptotic waves; one is an approximate Coulomb wave with p in the 
middle-range region, and the other a free (no-Coulomb) wave in the longest-range region. 
The constant Wronskians of this Hamiltonian can be calculated in either region. We can 
evaluate the Coulomb-modified nuclear phase shifts for the screened Coulomb problem, using 
the matching condition proposed by Vincent and Phatak for the sharply cut-off Coulomb 
problem. We apply this method first to an exactly solvable model of the aa scattering with 
the Ali-Bodmer potential and confirm that a complete solution is obtained with a finite p. 
The stability of nuclear phase shifts with respect to the change of p in some appropriate range 
is demonstrated in the aa resonating-group method (RGM) using the Minnesota three-range 
force. An application to the pd elastic scattering is also discussed. 
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§1. Introduction 



| In the momentum representation, incorporation of the long-range Coulomb force 

always poses problems. In particular, three-body scattering problems involving the 
Coulomb force are still under intensive investigations {DJDJSDJUiGD Here we mainly 
£S) | consider a much simpler problem of solving the Lippmann-Schwinger (LS) equations 

!y-^ ' in the momentum representation, in which the Coulomb force is included in the two- 

cluster resonating-group method (RGM). In this particular case, the longest range 
direct potential consists of a nuclear direct potential and the long-range Coulomb 
potential in the error function form when simple harmonic-oscillator shell-model wave 
functions are employed for clusters. We introduce a sharp cut-off radius p for the 
^ \ Coulomb force acting between constituent particles. We can solve the LS equations 

and obtain the T-matrix in the standard procedure. A problem is how to extract 
the correct nuclear phase shifts from this T-matrix or the phase shifts, including 
the effect of the screened Coulomb force. Here we propose a simple method, taking 
examples of the aa RGM and the proton-deuteron (pd) elastic scattering using the 
quark-model baryon baryon interaction. 

The standard procedure to solve the Coulomb problem, including the short-range 
nuclear potential and the long-range Coulomb force is well established as long as 
two-body problems are concerned. In the treatment in terms of the distorted waves, 
the relative- wave function tpi(r) between two clusters are solved numerically in the 
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configuration space, including the complete Coulomb force. The nuclear phase shift 
8f is then obtained from the asymptotic form of the relative wave function through 
the so-called matching condition 

W[Ge,ipe] p 

where Ft = Fi(k,r) and G( = Gi(k,r) stand for the regular and irregular Coulomb 
wave functions, respectively, and W[f, g] = f(k,r)(d/dr)g(k,r)—g(k,r)(d/dr)f(k,r) 
is the Wronskian for functions f(k, r) and g(k, r). The Wronskian values in Eq. (jl-ip 
are evaluated at the relative distance r = p, which should be taken large enough to 
avoid the effect of the nuclear force in the short-range region. Quite naturally, this 
standard procedure should be modified in the momentum representation, if we try to 
solve three-body problems like the pd scattering, and also the Lippmann-Schwinger 
RGM (LS-RGM) equations with the Coulomb interaction. In these applications, 
the basic ingredient is the T-matrix, which is usually formulated in the momentum 
space. The Born term of the T-matrix is already singular for the diagonal part of 
the initial and final momenta, qj = q{. 

A practical method to deal with the Coulomb force in the momentum represen- 
tation is to use the cut-off or screened Coulomb force. In the early work by Vincent 
and Phatak,® the Coulomb force in the -k^ + 16 scattering is assumed to be a 
sharply cut-off Coulomb force 

Up (r) = ^e(p-r) , (1-2) 

where n = a/hv is the Sommerfeld parameter and 8 is the step function. Since 
the relative wave function has a Coulomb-free asymptotic behavior, the asymptotic 
wave is composed of the nuclear plus cut-off Coulomb phase shift and the known 
Bessel and Neumann functions. This phase shift is obtained by solving the LS 
equation for the T-matrix in the momentum space. The nuclear phase shift 5f is 
then calculated from the matching condition of the asymptotic waves: 

W[F t ,u e ] p + tan S p e W[F e ,v e ) 
tan Of = - =o , (1-3) 

W[G e ,u e ] p + tan % W[GtM P 

with the sufficiently large p. 

A recent Coulomb treatment by Deltuva et al.^ 1 ®''^' for the pd scattering uses 
a screened Coulomb potential in the form of 

u f r ) = ^e-Wrf" (1-4) 
r 

with n ~ 4, and the "screening and renormalization procedure" , which is developed 
by Alt et a lpI3,[D,m ,113,0 The basic ingredient of this approach is the Taylor's 
theoremp^'GSli which implies that the phase shift of the screened Coulomb potential 
5^ requires the renormalization 5^ — > ag — C p (k) as p — > oo in the sense of 
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distribution, where 



C P (k) = ^ / w p (r) dr 




(1-5) 



is the diverging renormalization phase determined from the explicit form of the 
screened Coulomb potential co p (r). Since the relative wave functions with the screened 
Coulomb potential always suffer the renormalization of this phase factor, the com- 
plete pd scattering amplitude is achieved only when the limit p — > oo is reached 
in the two-potential formula for the scattering amplitude. In practice, this limit is 
taken numerically such that the well converged result is obtained. A problem of this 
procedure is that the error estimate of the finite p is not possible, and we usually need 
to take a very large p value, for which solving the Alt-Grassberger-Sandhas equation 
(AGS equation accurately is not easy because of the quasi-singular nature of the 
screened Coulomb potential. The convergence of the partial wave decomposition also 
becomes problematic, if the Coulomb singularity is so strong. 

A final goal of this study is to find an approximate but practical method to 
incorporate the Coulomb force to the pd elastic scattering, by using a reasonable 
magnitude of p. For this purpose, we incorporate the Vincent and Phatak approach 61 
to the "screening and renormalization procedure". In this paper, we first consider 
a simple potential model for the aa scattering and examine if this approach gives a 
reasonable accuracy of the phase shift, using the screened Coulomb potential. The 
stability of the nuclear phase shift with respect to the change of p in an appropriate 
range is examined by aa LS-RGM. An application to the pd elastic scattering is 
briefly discussed. 

In the next section, we discuss the sharply cut-off Coulomb problem, for which 
analytic derivation of the cut-off Coulomb wave functions is feasible. The definitions 
of the pure Coulomb wave functions used in this paper are gathered in Appendix 
A. A general procedure to calculate the nuclear phase shift from solutions of the LS 
equations for the two-cluster T-matrix is discussed in §3. In §4, a formulation for the 
screened Coulomb problem is given, by paying an attention to new features appearing 
in the screened Coulomb potential. An extension to deal with the pd elastic scattering 
in the present approach is given in §5. In §6, numerical performance is examined, 
first for an exactly solvable model in the case of Ali-Bodmer's phenomenological aa 
potential, secondly for the aa LS-RGM using the Minnesota three-range force, and 
finally for the pd elastic scattering using the quark-model baryon-baryon interaction 
fss2. In Appendix B, shift functions of various screening functions are evaluated. 
The screening function a p (R) for the pd scattering is derived in Appendix C. The 
last section is devoted to a summary and outlook. 

§2. Exact solutions of the sharply cut-off Coulomb problem 



In this section, we assume a sharply cut-off Coulomb potential Eq. (jl-2j) and con- 
sider the pure Coulomb problem as the limit of p — > oo. The regular solution <Pp(k, r) 
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corresponding to the Jost solution satisfies the following integral equation.® jUStjEDj 

1 f r , 2kr] , , , 

Pe( k > r ) = JITiM kr ) + J Goe(r,r';k)—6(p-r')ip p e (k,r')dr'. (2-1) 

Here, ui(kr) is the Riccati Bessel function and the Green function Qoi(r,r';k) is 
given by 

Goe(r, r'; k) = ^ [ue(kr) ve(kr') - ve(kr) ue(kr')] 8(r - r') , (2-2) 

with vi{kr) being the Riccati Neumann function. For r < p, ip p (k,r) is the same as 
the regular Coulomb function <pe(k, r) given in Eq. (|A-ip . which implies 

for r < p , (2-3) 

where F p {k) is the Jost function of the sharply cut-off Coulomb potential and F^{k) 
the Coulomb Jost function. If we use this in the integral equation for the regular 
solution ip%(k,r), 

r) = ~u e (kr) + {r\Gotu p i%) , (2-4) 

we obtain 

F p (k) 1 

ipe(k,r) = £ -j^u e (kr) + (r\Go£U} p ipt) for r < p , (2-5) 

with Gqi being the regular Green function. The Jost function Fp{k) of the sharply 
cut-off Coulomb potential is calculated from F p (k) = 1 + k e (ujg ^\uj p \ip p ) as^J 

F P(k) -■ - V {l + n)\ (-2ikp) m r{n + m-irj) 

&[) ~ %J] r{l+ 1 -in) ^ n\ ^ ml r(l + n + m + l) 

n=0 m=0 

xF(n + m — irj, £ + n + m + 1, 2ikp) . (2-6) 

In particular, 5"-wave Jost function is very simple: 

FS(k) = F{-i V ,l,2ikp) . (2-7) 

If we use the asymptotic form of the confluent hyper geometrical function F(a,j, z) 
at \z\ — > oo, we obtain 

«m (2k P )-'-Ff(k) = FHk) = , (2 . 8) 



resulting in 



lim ^§-(2kp)^ = 1 . (2- 



A Practical Method to Solve Cut-off Coulomb Problems 



5 



This relationship yields the limit of Eq. (|2-5p as 



Mk,r)= lim \ju e (kr)(2k P r + (r\GoeoJ p A)\ . (2-10) 

Furthermore, Eq. (|2-3p implies 

lim (2kp) iri ip^(k, r) = ip e (k, r) (2-11) 

p— >oo 

The asymptotic behavior of the cut-off Coulomb phase shift for p — > oo can be 
derived from the non-Coulomb version of Eq. (jl-ip . since Wronskians W[ui, ipe\ p and 
tpelp with the p — > oo limit are analytically calculated. The result is, of course, 

5p — > d£ — 77 log 2 as p — )• 00 , (2T2) 

with the ambiguity of integral multiples of 71- fiSl* ^ESJ 

The Jost solution for the sharply cut-off Coulomb potential is defined by the 
integral equation 

2kr] i 



f^(k,r)=J e + \kr)+ [ goe (r,r';k)^9(p-r')f^k,r')dr' , (2-13) 



where the Green function is 

g e{r, r'; k) = -i [u e (kr) ve(kr') - ve(kr) u e (kr')] 9(r' - r) . (2-14) 

The asymptotic behavior is given by 

fj?(k,r) = J+\kr) ~ e i(kr-(n/2)e) for r > p ( 2 . 15 ) 

For the Coulomb solutions, we cannot formulate the integral equation, since the 
asymptotic behavior is different from Eq. (|2T5p , However, for r < p, f p (k,r) can be 
written as a linear combination of two independent Coulomb Jost solutions, fi(k,r) 
and ft(k,r): 

j^k, r) = C p f e (k, r) + C£ft(k, r) for r<p. (2-16) 

The coefficients, Cf and C% are derived by evaluating the Wronskians W[fjt , /?] and 
W[/fc/f] at p -too. We find 



C ?=( 1 -^) (2M "' 



C p = (-Y+ 1 ^-(2kp)- iri e 2ikp as p -> 00 . (2-17) 
2kp 

Thus, if we use the symmetry Eq. (|A-5P for f£(k,r), we find 

f t P (k,r) = (l-Jj-\ (2kprf e (k,r) - ^L e ^(2kp)-^e 2ikp M-k,r) 

~ (2kp) ir 'f e (k,r) for r<p^oo. (2-18) 
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After all, we have obtained 

lim (2kp)-^fP(k, r) = fi(k, r) . (2-19) 

p— »oo 

Note that the renormalization phase is the complex conjugate of the one appearing in 
Eq. ()2-lip . This results in a basic property of the sharply cut-off Coulomb potential 
that the Coulomb Green function can be obtained as the p — > oo limit of the sharply 
cut-off Coulomb Green function. Namely, if we define 

G p e (r,r';k) = -r e (k,r < )fP(k,r > ) 

Gf(r,r';k) = -Mfc, *•<)£(*;, r>) , (2-20) 

then we find 

lim G p Ar, r';k) = G?(r, r'; k) for r, / < p oo . (2-21) 

This relationship is valid only when the Green functions are operated on the short- 
range potentials. 

One can derive the Coulomb scattering amplitude from the scattering amplitude 
for the sharply cut-off Coulomb potential. We use the formula for the short range 
force 

fp l i i i /P\ l i \rrP\ \ l$*mFg(k) 
ft = -p u e\u P m) = -js{ut\T£\u t ) = FP 

and calculate 

ft = lim (2k P yv f t (2kp) ir > ■ (2-23) 

p— »oo 

Equation (|2-8p yields for p — > oo 

1 , 0„„. .X 1 



aa ^-l)-»(«*-l)-Jf-»(<*rt"' 



Here, the last term is ^-independent and contributes only to 9 = if we add up over 
all the partial waves. Thus, we find 



/(6») = ^(2£ + 1)/ / P<(cob 9) 

oo 1 

= f{0) - £(2* + ^2ik ( {2kp)2lV ~ X ) Pe{c ° S 9) 

1=0 

= f c {6) for 0^0. (2-25) 
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Here, 

1 oo 



2ik 



V e ~2ir, log (sin |) A 1 + jg) 

2fc(sin|) 2 r(l-ir/) 



(2-26) 



is the standard Coulomb scattering amplitude. 

We should note that the renormalization phase {2kp) tr) appearing in the above 
equations is nothing but the Taylor's phase factor Eq. (|1 5|) . In fact, we can easily 
show for Eq. ([F2")l 



i r°° f p 1 

C(k) = — I io p (r) dr = V I - dr = V log {2k p) . (2-27) 

J 2k J 2k 

Then, the relationship in Eqs. (|2Tip and ()2T9p can be written as 

lim e iC "( fc )<(A:,r) = <ip t (k,r) , lim e -* P W f?(k,r) = f e (k,r) . (2-28) 

p— >0O p— >oo 

More basically, different parametrizations of Coulomb functions in Eq. (|A-9j) and a 
trivial relationship 

^ +) (fc,r) = ishn/^(fc,r) + /f f t (k,r) (2-29) 

derived from them are essential. Since many relations are also valid even for more 
general screened Coulomb functions introduced in §4, we reformulate the sharply 
cut-off Coulomb problem in more general form, using the parametrization of wave 
functions as 

tf(k,r) = ^e iS i F[(k,r) , 

rt(k,r) = ^\F p (k)\ F p (k,r) = real , 
f p (k, r) = e~^ [G p (k, r) + iF p (k, r)] , 

f p *(k,r) = e iS '[G p (k,r)-zF p (k,r)] . (2-30) 

For the sharply cut-off Coulomb force, the basic screened Coulomb wave functions 
satisfying 

lim F p (k,r) = F t (k,r) , lim G p Jk,r) = GAk,r) , (2-31) 

p— >oo p— >oo 

for a fixed r are analytically derived using the regular and irregular Coulomb wave 
functions, Fi(k,r) and Gi(k,r), and various Wronskians between these wave func- 
tions and the free wave functions. They are given by 



F p Ak,; 



ui(kr) cos 5p + vi(kr) sin 5 P „ r > p 



G p Jk,: 
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' jg^^r) + I (-W[Gi(k, r),u e (kr)] p sin % 
+W[G £ (k,r),v e (kr)] p cos 6 P ) F e (k,r) 
Vi(kr) cos 5p — ui(kr) sin b p t 



r < p (2-32) 
r > p . 



The screened Coulomb wave function tp p (k,r) also has a expression similar to Eq. 
(222]): 



(k, r) = -9m //(*, r) + // r) , 



(2-33) 



where /f = (l/k)e lS e sin 5?. Since f p (k,r) = Up(kr) = v^{kr) + iui(kr) for r > p 
the asymptotic form of Eq. (|2-33p is 



t/j p ^ + \k,r) = —ue(kr) + f p uj^' (kr) for r>p 



f p, .(+)/ 



We multiply Eq. (12^33]) by e^ P W and find 
< {+) (fc,r)e^ 



^m {//(A;, r)e-^} 



+ 



a <C'(fc) fPJC{k) 



f p (k,r)e-^ 



where we have set 



fP = 

Jr) — 



2ik 



(2-34) 



(2-35) 



Here, we take the limit p — > oo and use Eq. (|2-28p . If we compare the resultant 
expression with Eq. ()2-29j) . we find the correspondence 



JC(k) fPJC(k) 



(2-37) 



In fact, Eq. (|2-37p diverges, but if we add up over all the partial waves, the second 
term of Eq. ()2-37|) does not contribute except for 6 = because f(] is I- independent. 
Thus, the scattering amplitude of the sharply cut-off Coulomb force 



ne) = Y / (2£ + l)f p P l (9) , 



satisfies 



lim e i(P{k) f p (e)e iC:P{k) = f c (0) for 0^0. 



(2-38) 



(2-39) 



§3. Two-body Coulomb problem 



In this section, we consider a scattering problem for a two-body Coulomb system 
consisting of a short-range local potential v(r) with the interaction range a and the 
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Coulomb force co c (r) = 2kt]/r. The Schrodinger equation in the configuration space 
reads 



e(£ + i) 



n 2A77 2 

ufrj h fc 

r 



with the boundary condition 



(+), 



r) ~ - c 3mf((k,r) + f e f e (k,r) 



M +) (r) 



r — )• 00 



, 



(3-1) 



(3-2) 



Here, r) is the Coulomb Jost solution in Eq. (|A-3p and the partial-wave scat- 
tering amplitude fg is expressed as 



with 



1 



1) 



(3-3) 



using the nuclear phase shift 5g . In the usual approach, 5f is calculated from the 
real regular function J-^{k,r) for the Schrodinger equation Eq. (|3Tp . which satisfies 
the relationship 



w\ + \r) = ~e< a * +5 ^T e (k,r) 

Ki 



(34) 



The asymptotic wave of 7^(fc,r) is expressed as 

Ti(k, r) ~ F/(fc, r) cos tff + G/(fc, r) sin 5f 

~ sin (At — 77 log 2 At — (tt/2)£ + an + <5j 

The nuclear phase shift 5^ is then calculated from Eq. (jl-ip by assigning J-g{k, r) to 
and taking p > a large enough. 



5f) 



00 



(3-5) 



Similar equations are also valid for the sharply cut-off Coulomb force 



Eq. (|l-2p . Namely, the Schrodinger equation for this system, 
2 



has the asymptotic wave 



+ D 



f (r) — 0Jp{r) + fc 



8£ w (r) = -9ta//(fc,r) + ##(fc,r) 



fltf (+) (r) = 



for r > a 



u) n \r) m 



(3-6) 



T uAkr) + fpbji + \k,r) for r>p 



(3-7) 



where f^(k,r) is the Jost solution for w p (r). The scattering amplitude /? in Eq. 
(|3-7p this time is parametrized as 



1 



Jiff 



and ff N 



1) 



2i8 p . f pN 



with ^ = $ + 5f , 



(3-8) 
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where /f = (l/2ik)(e 2lS i — 1) is the scattering amplitude for oj p (r). Furthermore, we 
have the relationship 



Fp(k, r) = Ui(kr) cos 5« + vi(kr) sin 5« (r > p) . (3'9) 



Note that the second equation of Eq. (|3-9p is exact for the sharply cut-off Coulomb 
force. We multiply Eq. (|3-7f) by the phase factor e^ P ^' and take the limit p —> oo. 
Then, a procedure similar to Eq. (|2-35p leads to the correspondence 



lim 



iC(fc) fPJC{k) 



p— ¥00 L 

lim ^ p(+) (A;,r)e iCP(fc) = <P ( f +) (k,r) for r>a. (3-10) 

p— >oo 



From Eq. (|3-8p . the nuclear phase shift 5f is obtained through 



lim 5f = lim (5£-<tf) . (3-11) 

p— >oo * p^oo * * ' 

The sharply cut-off Coulomb phase shift 5 P is calculated from 

Since <5£ is obtained by solving the potential problem for v(r) + uj p (r), Eqs. (|3T1|) 
and (|3T2p gives a solution for the two-body Coulomb problem in the momentum 
representation, using the sharply cut-off Coulomb force. 

Another method to derive the nuclear phase shift in the momentum represen- 
tation is to use the two-potential formula for the T-matrix. For the short-range 
potential v and the sharply cut-off (or screened) Coulomb potential ui p , we solve the 
T-matrix equation 

T p = (v + lu p ) + (v + oj p )G T p , (3-13) 

where Go = (z — /io) _1 with z = E + ie is the free Green function with the energy 
E. We assume the energy factor (h 2 /2fi) = 1 and set E = k 2 . Furthermore, the 
partial wave decomposition is implicitly assumed and the orbital angular momentum 
£ is omitted for typological simplicity. The kinetic energy operator ho is, therefore, 
ho = (d/dr) 2 — £(l + l)/r 2 . The two-potential formula for T p is given by 

T p = t Up + (1 + t Up G %(l + GoW 



tui p = w p + UpGotujp = uip + u} p G Up ujp 

t Up = v + vG Up t Up = v + vG p v , (3-14) 

where G Wp = (z — ho — and G p = {z — ho — v — oo p )~ l ■ To derive the scattering 

amplitude, we sandwich T p with the plane wave (with the wave number k) 

\<t>) = \W) , (3-15) 
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and define 

= \(f>) + G cu p \r {+) ) ■ (3-i6) 

Then, by using 

= (l + Go^ p M , (r { - ) \ = (<f>\(l + t l , p G ), (3-17) 

we find 

<0|r>|0> = <«^» + <^ H M< (+) > . (3-18) 

In Eas. (I3~T7)1 and (l3~T8]l . is defined by ^~\k,r) = (A;, r))*. This 

equation is essentially equivalent to the T-matrix in the distorted-wave Born ap- 
proximation (DWBA). In fact, if we set 

t^\V {+) ) = v(l + G Wp t Wp W (+) ) = v\^) , (3-19) 

the on-shell T-matrix is expressed as 

{4>\T P \<P) = (0|w p |V p(+) ) + (V pH M^ p(+) > • (3-20) 

The LS equation for the total wave function 

|SW(+)) = \^p(+)) + G Up v\9^) (3-21) 

is equivalent to Eqs. (|3-6|) and (|3-7|) . 

Equation (|3T8p gives a starting point for the "screening and renormalization 
procedure" . Namely, if we sandwich Eq. (|3T8p with the renormalization phase e l ^ P 
with (, p = C p (k), and take the limit p — > oo, we find 

UlTU) = lim e^'MT^e*' 

p— >oo 

= lim e^^lt^e^ + lim e'WV^I&^+V^ 

p— >oo ^ ' p— >oo 

= (<^ c |0) + lim(^(-)|t£|^(+)) . (3-22) 

p— >oo 

Here, 

I^W) = lim |V£ {±) }e ±<p (3-23) 

p— >oo 

are the pure Coulomb wave functions. The first term in Eq. (|3-22p is separated into 
the partial-wave Coulomb amplitude ff and ^-independent term from the discussion 
of the preceding section. When all the partial-wave contributions are added up, the 
first term becomes the pure Coulomb amplitude. Actually, the relationship between 
the scattering amplitude and the on-shell T-matrix yields 

WcM=~ 3 f/ c m=~^ c m (3-24) 
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with \qj\ = |qj = k. In the second term of Eq. (|3-22p . p — > oo limit can be 
taken, since the nuclear potential v is short-ranged. We define t by the solution of 

t = v + vG c t, where G c = (z—ho — Irjk/r) -1 is the Coulomb Green function. Thus, 
we find 

lim (^M^) = . (3-25) 

p— >oo 

To derive this matrix element, we introduce the total wave function l^^) through 

tj^(+)) = V {1 + G c t)\^ +) ) = v\^) , (3-26) 

which satisfies the LS equation 

|#(+)) = (l + G c t)\ipM) = |^ {+) ) + G c v\&^) , (3-27) 

and the Schrodinger equation in Eqs. (|3Tp and (|3-2p . We should note that Eq. ()3-27p 

has a solution, since t> is short-ranged. Here, we introduce a decomposition of the 
partial-wave Green function 

G C (r,r';k) = -Mk,r < )f e (k,r > ) , 

G c = G c - \ft){ip i f ) \ with G c ^0 as r -)• oo . (3-28) 

Then, we find the asymptotic behavior 

= |^ (+) ) - |/)(V> H M^ (+) > + G c ^ (+) ) 

~ - l/)<^ (_) l*l^ (+) > as r^oo. (3-29) 

If we use the Wronskians of the Coulomb wave functions 

W[F e ,G e ]=-k, 
W[F e ,f e ] = -ke- i<T ' , 

W[4 + \f e ] = -l , (3-30) 
derived from Eq. ()A-9p . we obtain 

W[hM +) ]r^oc = 1 (3-31) 

and 

Wl), = Ita ^If (*.r).f;»(*.r)l (3 . 32) 

If we further parametrize 

^-)\t\^) t = -e 2 ^^ (e 2 ^ - l) , (3-33) 
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Eq. (|3-32p is equivalent to 



_ _ W[W,r) ^<+>(fc,r)] 
1 r-4oo W[G t (k,r),VM(k,r)] y ' 



Eventually, we find 



= /<?(0) + £(2^ + l)e 2 ^ /f P,(cos 



£=0 

with /f = ^ (>f _ i) . (3-35) 

For practical calculations in the momentum space, it is much easier to start 
with the sharply cut-off Coulomb force from the very beginning. We multiply the 
LS equation in Eq. (|3-2ip by the renormalization phase e l< ~- p and take the limit p — > oo. 
Then, by using Eqs. (|2-11|) and (|2-2ip . we obtain 

lim |^(+)) e <" = + G c v lim |t^(+)) e iC * . (3.36) 

p— >oo p— >oo 

If we compare this with Eq. (|3-27|) . we find 

lim \WrP(+)\jC P = |^(+)) _ (3.37) 

p— >oo 

If we further use Eq. QIES?D in tan df of Eq. (EM . we find 

tan „ _ _ ]im WW. r), »?+'(*■,)] 
"~»'[G,(fc,r),< +, (k,r)] 

for sufficiently large p. On the other hand, the asymptotic behavior of the wave 
functions for the short-range force yields 

wf + \k,r) = -e^ {ui(kr) cos S p e + v £ (kr) sin ^} 

for p < r — > 00 , (3-39) 

for sufficiently large p. Thus, if we calculate Wronskians in Eq. (|3-38|) at r = p, we 
obtain 

N W[F e (k,r ),tf {+) (k,r)] r=p 
tan di = 



W[G e (k,r),^ (+) (k,r)] r=p 
W[F e ,u e ] p + t!m6 p e W[Fe,ve] f 
'w[Gi,u e } + tan 6% W[G e ,v t ] 



(3-40) 



which is nothing but Eq. (|l-3p . After all, if 5 P is calculated in the momentum 
representation, the nuclear phase shift 5f is obtained through Eq. (|3-40j) . The 
scattering amplitude fe(0) is calculated from Eq. (|3-35|) . using 8f. 
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§4. The screened Coulomb case 

In this section, we will extend the preceding discussion for the sharply cut-off 
Coulomb force to a more general screened Coulomb force, which is formulated as 

Up (r) = ^%(r) , (4-1) 

according to TaylorP^ 1 Here, the screening function a p {r) with 1 > a p (r) > is a 
monotonically decreasing function of r, satisfying 

1) with p fixed, a p (r) decreases to zero, faster than 0(r~ e ~ 2 ) (e > 0), as r ap- 
proaches to oo, 

2) with r fixed, a p {r) appraoches to 1 as p approaches to oo, 

3) around r ~ p, there exist sufficiently wide regions in which a p {r) ~ 1 and ~ 0. 

In 3) above, we added "an almost sharply cut-off condition" in addition to 
the original conditions 1) and 2) in Ref. ll5|) . This condition is required if we wish 
to develop almost parallel discussion to the sharply cut-off Coulomb case, as seen 
below. Note that the sharply cut-off Coulomb case is included in the above category 
by taking a p (r) = 9{p — r). 

The necessity to relax the sharply cut-off condition is as follows. First, in the 
LS-RGM formalism, the longest-range direct Coulomb potential becomes a screened 
Coulomb force as explicitly shown in §6.1 and §6.2. If the cluster wave functions are 
assumed to be standard harmonic-oscillator shell-model wave functions, the cut-off 
function ot p (r) is usually expressed by the error function. Secondly, in the application 
to the pd elastic scattering, the asymptotic Hamiltonian involves a screened Coulomb 
force which is obtained from the pp Coulomb force by the folding procedure using 
a realistic deuteron wave function. In Ref.[7j), the same pp screened Coulomb force 
is used for the pd screened Coulomb force, but using a more realistic pd Coulomb 
potential is certainly desirable to avoid unnecessary extra distortion of the deuteron 
in the asymptotic region by the Coulomb force. In any case, the screening func- 
tion CKp(r) should be chosen most appropriately for each problem, since "in practice 
Coulomb potentials are always screened" as stated in Ref- fTo"]) . 

For the screened Coulomb force in Eq. (|4-ip . the parametrization of screened 
Coulomb wave functions in Eq. (|2-30p is employed in the following, but the explicit 
solutions of F{?(k,r) and G^(k,r) like in Eq. (|2-32|) are no longer available. In order 
to extend Eq. (|2-32p to the screened Coulomb case, we first examine the behavior of 
the screened Coulomb wave functions around the origin r ~ 0. For the pure Coulomb 
solutions, Fi(k,r) and G^(/c,r), we can easily show that 

*>( fc ' r ) ~ T^JTs-Mkr) , G t {k,r)~\F t (k)\v t {kr) as r^O, (4-2) 

by using the explicit expression of the Coulomb Jost solution fi(k,r) in Eq. (|A-3|) 
and the parametrization Eq. (|A-9p . The corresponding expressions for the screened 
Coulomb wave functions are 

... . 1 
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G p e {k,r)~\F[{k)\v e (kr) + — ^— A p (r)u e (kr) as r^O, (4-3) 

where an extra term including A p {r) appears in the irregular solution G p (k,r). The 
real function A p (r) is given by 

A p (r) = -\F p (k)\±-W[G p (k,r),v e (kr)] 

{log r £ = 

for as r -> , (4-4) 

and diverges as r — >• 0. These results are derived by applying the Calogero's variable 
phase method^ to the regular solution ip p (k,r) and the Jost solution f p (k,r). 

For practical applications, we use the "almost sharply cut-off condition" 3) and 
assume a screening function satisfying 

3)' 

fl for r < p — b = R m 
Q p (r) = ^ , (4-5) 

[0 for r > p + b = R out 

with a sufficiently large p 3> b. A new parameter b is introduced to make a smooth 
transition for the Coulomb force to disappear. To make the pure Coulomb region 
available, i?j n S> a should also be taken large enough, compared with the range a 
of the short-range nuclear force. By this assumption, we can extend the discussion 
in the sharply cut-off Coulomb case, although some modifications are necessary as 
seen below. First we apply Calogero's variable phase method to the regular function 
(p p (k,r). This solution and the pure Coulomb wave function (pi(k,r) both satisfy 
the same integral equation Eq. (|2-ip with 6{p — r) — > a p (r) or 1 for r < R m , yielding 

(p p (k,r) = (pe(k,r) for r < R in . (4-6) 

If we use the standard relationship 

Mk,r) = ^ I \F e (k)\F e (k,r) , (4-7) 



(see Eqs. (fA^ and $FW\\ ). Eq. implies 

\Fi(k)\ 



F p {k,r) = S,,J, ifr(fc,r) for r < R in . (4-8) 



Here, we can prove 



S = 1 • (4-9) 
P^oo \F p (k)\ v ' 



Similarly, the screened Coulomb phase shift 5 P = 5 p (k) is proved to have the Coulomb 
limit 

lim (5 P + ( p ) = a e , (4-10) 

p— 5>00 ' 
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where C, p = ( p (k) is given by Eq. fjl - 5[) for 0J p (r) in Eq. ()4-ip . For the irregular 
solution G p (k,r), the local phase approach does not work. In this case, we have 
admixture of the regular solution Fg(k,r) for r < i?j n , which is related to A p (r) in 
Eq. (|4-4p . Summarizing the above discussion, the explicit results of Eq. ()2-32[) in the 
sharply cut-off Coulomb case should be modified to 

F p (k,r) = ^F e (k,r) , 



G p e (k,r) 
F p (k,r) 
G p e (k,r) 



a p G e (k, r) + A p F e (k, r) for r < R h 
ui(kr) cos 5 P + Vi(kr) sin S p , 
vi(kr) cos 5 P p — ue(kr) sin 5 P „ 



for r > R, 



out ; 



where 



and 



\F e p (k)\ 

\Fm\ 



lim a p f 

p— >oo 



(4-11) 
(4-12) 



as p 



with 



C = C p (k) 



1 

2k 



oo 



u! p (r) dr . 



(4-13) 



We note that, for the pure Coulomb problem, the renormalization of the screened 
Coulomb wave functions and the scattering amplitude is possible. In particular, 
Eqs. (|2 33|) - (|2-39p are all valid owing to Eq. (|4-13p . However, the renormalization of 
the irregular solutions like in Eq. (|2-19|) needs a modification, since in general A p ^ 
in Eq. (|4-lip . For example, the relationship in Eq. (|2-28p should be modified as 

lim ~ e ^+C) F p (k-r ' 



lim 4>f + \k,r)e i(P 

p— >oo 



,r) = -e 
k 



i a, 



p— S>00 k 

(for regular Coulomb wave function) , 
lim f p (k,r)e-^ P = lim e^K") \G p (k,r) + iF p (k,r 



F e (k,r) 



V^ +) (£:,r) 



p— >oo 



GAk,r) + lim A p F t (k,r) + iF e (k,r) 

p— >oo 



h(k,r)+A e e-^ F e (k,r) 

1 



lim 

p— >oo _ 



2ik 



(e*> - 1) = fj 



(4-14) 



for r < R[ n . Here, we have assumed that Ai = lirrip^oo A p exists for simplicity. By 
the same token, the p — > oo limit in Eq. (|2-35p becomes 

^ e + \k,r) = {f e (k,r) + Ate'™* F e (k,r)} 
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. /'i lor /• < liin ■ (4-15) 

Here, because Ag is real, the contribution from the terms proportional to Ag vanishes 



as 



-Ag sin at F t 
k 



-l(T e 



Fg(k,r) = 



(4-16) 



resulting in Eq. (|2-29j) again. It is important to note that this renormalization 
is possible only for the regular solution of the pure Coulomb problem. Once the 
nuclear potential is introduced, we need further renormalization for the magnitude 
of the wave function related to A p , since the derivation of the regular solution also 
requires irregular solution of the screened Coulomb problem. 

In order to make the similarity to the sharply cut-off Coulomb case more trans- 
parent, it is convenient to introduce a modified set of screened Coulomb wave func- 
tions by 

F p (k,r)=a p F p (k,r) , 

G p (k, r) = - p G p (k, r) - A p F p (k, r) . (4-17) 
a e 

For r < Ufa, these are the pure Coulomb wave functions: 

F £ p (k,r) = F e {k,r) , G p (k,r) = G e (k,r) for r < R in . (4-18) 
However, for r > -Rout, Eq. (|4-lip leads to 

F p {k, r) = a p \ug(kr) cos 5 P + ve(kr) sin 5 P ~\ , 

~ 1 

G p (k, r) = — [vi(kr) cos 5 P — ue(kr) sin 8 P ] 
a i 

—A p [uf(k, r) cos 5 P + ve(k, r) sin S p ] for r > R out . (4-19) 

The Wronskians of these wave functions with the free scattering solutions in the 
asymptotic region are given by 

1. 
k 
1, 



-W 
-W 



k 



1. 

k' 



-w 



Fj?(k,r),u e (kr) 
F p (k,r),v e (kr) 
G p E (k,r),ut(kr) 

G p Jk,r),v e (kr) 



a p sin S p , 



-a i cos 



Of , 



— cos 5 P — A p sin 5 P , 

do 



sin 8 P , + A p cos 5 P for r > R out . (4-20) 



Let us assume a <C R{ n and consider the regular solution of the Schrbdinger equation 
for v(r) + (2krj/r)a„(r)'. 



^ + \ r ) = F p (k,r) cos 8f + G p (k,r) sin 5? 



for r > a 



(4-21) 
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In the a < r < Ri n region, 8f becomes the nuclear phase shift owing to Eq. (|4T8p , 
This can be calculated from 



tan 5, 



N 



w 



F p {k,r),^ 



p(+), 



W 



for r > a . 



(4-22) 



The Wronskians in Eq. (|4-22p can be calculated at any points r > a, since F p (k,r), 
G p (k,r) and !?f(r) are all solutions of the Schrodinger equation for the screened 
Coulomb potential. In particular, the asymptotic behavior 



p(+) 



(r) = B ui(kr) cos 8 P + vi(kr) sin 8g 



for r > R t 



out 



without the Coulomb force, yields a connection condition 



Ti- 



tan 8 



N 



F p (k,r),u e (kr) 



+ tan 5 P P W 



Fj?(k,r),vt(kr) 



w 



G p Jk,r),U£(kr) 



Rou 



+ tan 5 P P W 



G p Jk,r),v e (kr) 



(4-23) 



(4-24) 



which is an extension of Eq. (|3-40p in the sharply cut-off Coulomb case. The phase 
shift 8 P is calculated from the standard procedure to solve T-matrix of v(r) + 
(2krj/r)a p (r) in the momentum representation. 

To the contrary, we can also recover the asymptotic behavior of !Pf^(r) in 
Eq. (|4-23p . starting from Eq. (|4-2ip and Eq. (|4-24p . If we use the expressions of 
Wronskians in Eq. (|4-2U|) . the connection condition Eq. (|4-24|) yields 



tan 5, 



N 



a p sin ( 8 



S p 



COS S p 5 p )+A p sin (8, - 8 p e 



(4-25) 



We write this as 



8111 Og 



cos 6 l = 



B p 



sin [So — 8 P 



1 



s cos [ft -8", )+A p sin [S'l - S p 



B p 



— cos (d^ — 8 p ^j + A p sin [s P t 



8 P 



+ 



a„ sin 5 



TP 



5 P 



(4-26) 



If we use this in Eq. (|4-21|) for r > R on t, the asymptotic behavior of Fp(k,r) and 
G p e (k,r) in Eq. (ICTjl yields Eq. (l¥23l) with B = l/B p . In particular, if v(r) = 0, 



8 P = Sp in Eq. f|4- 25j) yields the correct results Sf = 0. 



In fact, Sf in Eq. g2B is p-dependent: 5 



N 



limit Sf = lim^oo 8% . Furthermore, the present assumption that a p (r) = 1 or 



and we need to take the 
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except for the interval [.Rim -Rout] = [p — b,p + b], is just an approximation. We 
have to examine the accuracy of this approximation for the finite p on the case- 
by-case basis. In practical calculations, we solve F?(k,r) and G p Ak,r) from R in to 
-^out; by taking the starting values of the pure Coulomb wave functions Fi(k, Ri n ) 
and Gi(k, Ri n ). The Wronskians needed in Eq. (|4-24p are calculated numerically. In 
the sharply cut-off Coulomb case with 6 = and i?i n = R out = p, this process is 
unnecessary, and reduced to Eq. (|3-40p . 

The extra term proportional to A p in Eq. (|4-17j) also affects the relationship 
of the Green function in Eq. (|2-21|) . To find a new relationship for the screened 
Coulomb force, we solve Eq. (|4-17j) inversely and express F p (k,r) and G p (k,r) as 



F p (k,r) = ^F p (k,: 



G p (k, 



'G p (k,r) + A p F p (k,; 



(4-27) 



Then the Green function of the screened Coulomb force in Eq. (|2-20p is expressed 
for a fixed £ as 



1 1 



G?=G? J -~- p -A p \F p )(F p l , 



with 



G p (r,r';k) = -l^F p (k,r < : 

K (X is 



a?G?(fe,r>)+*^F/(fc,r>; 



For r, r' < i?; n , the p — > oo limit of Eq. (|4-29p yields 



(4-28) 



(4-29) 



lim G p (r,r';k) = Gf(r,r';k) - T A e \F e )(F e 
We keep the finite p and write Eq. (|3-2ip as 



for r, r < R n 



|#p( 



+)\ 



x 



+ G p v |^+)) - ~ - p A p t \F p )(F p \v\W p ^) 

rC Qj£ 

I e iS "^ \F P ) \l - e- iS iA p AF p \v\^P^)] + G p u v 



+)\ 



(4-30) 



(4-31) 



Here, we define 



Then, we find 



\\pp(+)) = |#p(+)\ 



iS tA p {F p 



|^(+)) = I J_ \ F P) + GPjV 



(4-32) 



(4-33) 



Here, we multiply Eq. (|4- 33|) by e^ P and take a limit p — > oo with r G [a, R- m ] fixed 



The first term of the right-hand side of Eq. (f¥33"j) is {l/k)e iue \F t ) = |^ +) ). In I he 
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second term, we further use the decomposition of the Green function 
1 1 



G p = 3L 



G p + i(l/a p )F p )(Ff\ 



g^(r,r';k) = -± {F p (k,r) G p (k,r') - G p (k,r) F p (k,r')} 9(r' - r) , (4-34) 



and find 



o 

1 1 

a' 



dr' GP(r,r';k)v(r , )W < - + \r') e^" = I dr' $f(r, r'; k)v{r')^ + \r') ^" 



k » p 



a p G p (k,r)+i-F p (k,r) 

Qjn 



dr' F p {k,r') v{r') ^ + \r') . (4-35) 



Here, the first integral in the right-hand side vanishes since v(r') = for r' > r > a. 
In the second integral, the range of v{r) makes r' < a only, so that we can safely 
replace F p (k,r') by Fg(k,r'). Thus, we find 

G^v |^ (+) ) e iC -" ~ -\\ [G t + iF e ]}(F e \v\V p ( +) } e K " 



le^\f e )(F e \v\V p M) e*' 
-|//>(^ (-) |w|^ p(+) > e K " for p 



oo , 



(4-36) 



and 



(+)) e *C" = ~ \ft) lim l^-hv\^ + h e K " 

p— >oo ' p— s>oo 



lim \<F p{ - 



for a < r < R in . (4-37) 



If we compare Eq. (|4-37p with the asymptotic form in the exact Coulomb case in Eq. 
(HE21, we find 



lim |^( + )) e i(P = \V { 



+)\ 



p— >oo 



(4-38) 



In the matrix element of Eq. (|4-32p . we can also replace F p by F%, since v(r) is 
short-ranged. By solving Eq. (|4-32[) inversely, we can show that 



\&ti+)) = |«p/>(+)) \i + e -iS p eA P (F e \v\V p(+) ) 



(4-39) 



If we multiply Eq. (|4-39p by e l< * P and take the limit p — > oo, 5 P — > <T£ — £ p yields 



lim |^+))e^ P = |<f (+) ) l + e- iae A e (F e \v\& {+) ) 

p— >oo L 



(4-40) 



and 



|<fW) = lim e iC " 1 - e- iae A p (FAvW p ^) e iC -" 

p— >00 



(4-41) 
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This expression implies Eq. (|3-37p is no longer valid for the screened Coulomb force, 
and we need an extra normalization factor [l — e~ tae A p {Fi\v\ } I ,p<y+ ' > ) e*^H 

Finally, we will show that another type of the connection condition, equivalent to 
Eq. (|4-24p . is also obtained by considering two types of asymptotic forms of \^F P ^). 
First, the asymptotic form of \\pP( + )) for R m > r — > oo is from Eqs. (|4-33p and (I4-35|) 



|^(+)) „ 1 J_ \ F P) 



1 l 

k a p D 



a p G p {k,r)+i^F p (k,r) 



{Fi\v\& p( - +) ) for r < R h 



(4-42) 



The Wronskians at r — > oo with r < R\ n are given by 
W 



(F e \v\^) = e i5 c -1 (F t \v\wf + ^ » p »~ iS ° 



a e e 



W 



1 



1 



K) 



- p {l-i{F,\v\^)- p e-^\ . (4-43) 



Thus, if we define K p „ by 



K p ±(F t \v\¥*+'>) 



p —i 



ar c e 



l-i(F e \v\<F p( - +) )-^e- 



we obtain 



K p W 



k W 



(4-44) 



(4-45) 



Here, we note that all the wave functions with tilde satisfy the Schrodinger equation 
for the screened Coulomb potential for r > o, so that we can evaluate Wronskians 
at any points r > a. If we take the limit p — > oo in Eq. (|4-44p . Eq. f|4-38|> yields 



lim l(F e \v\* p{+) He- i5 ' = l(F e \v\^ {+) )e- iat = Uf^^)^ 

p— >-oo fx, Ki Ki 



^{F t \t e \Fi) = e- 2i(Tl l^ i \iA l - i + ] 



Thus, if we define lim,^*, K p = K? Eq. flUi becomes 



>,Ri7j„/,(+)\ 



If we further parametrize 



we find K, 



N 



(4~ ] \te\4 +) ) 
-k cot . 



-c- ^ — !— ( e 2l& ^ - I 
2ik 



2ia e 



(4-46) 



(4-47) 



(4-48) 
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On the other hand, in the region r > R Q ut, the Coulomb-free asymptotic wave 
gives 

= l\ ue ) - |4 | - ) }(&|27|&> for r > # out , (4-49) 

where the T-matrix Tj? is defined in Eq. (|3- 13[) . If we write Eq. (|4-49j) by the K- 
matrix defined by 

RP (<l>e\Tf\<t>t) = l-ik (falTflfa) , (4-50) 

it is expresses as 

|SW*+>) = [\ut)K p t - \v e )k] l(<Pe\T e p \<P e ) for r > R out . (4-51) 

We can use this to calculate the Wronskians in Eq. (|4-45|) at r = R Q ut, since the 
difference between !^ + ^ and $^ + ^ is just a normalization. From these processes, 
we eventually obtain 

K$ \w[F p t MR^K - W\FlM*J\ = k [wfilMR^K* - kW[G p e ,v e ] Rout } , 

(4-52) 

which is equivalent to Eq. (|4-24p since K p = — k cot 5 P and K% = — k cot 5% . 

Summarizing this section, we first calculate (<j>e\Tt \<f>i) from the LS equation in 
the momentum representation, calculate K p by Eq. (|4-50[> . transform to K p by Eq. 
(|4-52p . and take the limit lim^oo K p = Kf . Then, the nuclear phase shift 5f is 
obtained from KF = — k cot8¥ . 

§5. Application to the pd scattering 

Application of the present formalism to the pd scattering is not straightforward 
because of several reasons. First, the asymptotic pd Coulomb potential suffers the 
strong distortion effect of the deuteron due to the long-range nature of the Coulomb 
force. In the strict three-body treatment of the nd scattering by the AGS equation, 
the distortion effect of the deuteron is fully taken into account, but only for the 
short-range force. Even if we neglect the Coulomb distortion effect by using the 
screened Coulomb force, the quasi-singular nature of this interaction causes the dif- 
ficulty that the treatment by the standard AGS equation eventually breaks down at 
the limit of p — > oo. To avoid this, a new formulation by the Coulomb-modified AGS 
equation was devised. However, very singular behavior of the screened Coulomb 
wave functions in the momentum representation makes it difficult to solve this equa- 
tion numerically. Another difficulty lies in the partial-wave expansion of the AGS 
equation. Even in the two-body Coulomb problem, the partial-wave expansion of 
the Coulomb amplitude does not converge in the usual sense, but converges only as 
the distribution. It is therefore attempted to formulate the AGS equation based on 
the three-dimensional description of the two-body t-matrixJ^®® 1 ®'® The isospin 
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symmetry breaking by the T = 3/2 component should also be taken into account, 
since the Coulomb force admixes the different isospins. Here, we extend the "screen- 
ing and renormalization technique" to incorporate the present approach and try to 
find a practical method to deal with the pd elastic scattering even in an approximate 
way. 

Let u p (r; 1,2) be a screened Coulomb force acting between two nucleons 1 and 

2: 

w'(r; 1, 2) = —a p {r) . (5-1) 

Here, r is the relative coordinate between the two nucleons. We use a set of Jacobi 
coordinates of particles (l-2)+3 as the standard one and denote it by 7 = 3. Another 
relative coordinate is denoted by R in this section. Then, the screened Coulomb 
potential Eq. (|5-ip is expressed as ojP with 7 = 3. In the following, we formulate the 
Coulomb-modified AGS equation in the isospin representation. The three-particle 
symmetric three-body screened Coulomb potential oJq = uja is given in the 
isospin basis as^ 1 

u p c = ojP + WP + WP for v 7 . (5-2) 

Here, W§ denotes the screened Coulomb potential between the nucleon 7 and the 
residual NN pair, and is a function of the Jacobi coordinate Ry between them. 
Furthermore, the three-body potential W7, which is usually called the polarization 
potentialpS is defined by 

w$ = E - 5 ^ w p) = E - w " ■ ( 5 - 3 ) 

p 

It should be noted that, for ppn system, either of ojP or WSj + WJy in Eq. (|5-2p is only 
non-zero. 

The two-potential formula for the three-body system is derived for the solutions 
of the Coulomb-modified AGS equation.i 1 ' First, the three-body transition operator 
Up a for the usual AGS equation is defined through, 

G p = Sp, a g p a + 9 p pUlX , (5-4) 

where the full resolvent Q p and the channel resolvent g p a are defined by 

9 P = ^z-H -J2v*-u£j , g p a = (z-H -v a -u } P a )- 1 , (5-5) 

with v a being the short-range nuclear potential and z = E + Ed + iO composed of 
the incident energy E and the deuteron energy e^. The three-body kinetic-energy 
operator is expressed as Hq = /io 7 + /107 for an arbitrary set of Jacobi coordinate 7. 
The transition operator U P n a satisfies the AGS equation 

^ = + E tP c °0 U',a » (5-6) 
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where Gq 1 = (z — -ffo) -1 is the free resolvent and basic two-nucleon T-matrix t p is 
generated by solving the LS equation for v a + ui p . Namely, 

jP = (y + u p ) + + UJ P ) G t p . (5-7) 

In Eq. (|5-6|) and below, we use the usual convention <5^ Q = 1 — 5p, a . The full resolvent 
Q p can also be decomposed as 

G p = Sp,aG p a + G p U P !a G p a , (5-8) 

using another resolvent Ga defined by 

G p a = (z-H -v a -u; p -W p )- 1 , (5-9) 

The operator Uo a satisfies the Coulomb-modified AGS equation^ 

K* = h,<* ((G p a y l + v a ) + 6p ja W p + Y, + W) G p U^ a . (5-10) 

(7 

From the relationship between g p and Ga, the operator U£ a is related to Up a 
through 

V% a = Sp,«T p + (1 + T p g p )U P ja (l + g p a T p ) , (5-11) 

where the screened Coulomb T-matrix T„ for the pd scattering is obtained from W£ 
through 

T p = W p + W p g p T p = W p + W p G p a W p . (5-12) 

Equation ()5- 1 1 j) is the two-potential formula for the three-body system. The Coulomb- 
distorted asymptotic wave function is defined by \ip$ ) = (1 + gaTa)\4> a ) from the 
channel wave function \<p a ) = IQoq'V'o)- Prom this definition and Eq. ()5T2p . we 
obtain 

\r a {+) ) = \<Pa)+g p a w p \r a {+) ) . (5-i3) 

We define \tpa ) as the complex conjugate of {ipa^} and find 

(M U U^ = 5 fi,M T M + {rt ] \UiM +) ) • (5-14) 
We can separate the deuteron part in Eq. ()5T3p and we obtain 

K (+) > = lx£ (+) ,^> , 

\x p a {+) ) = \Qoa) + (E a + i0- hoar'WMx^) , (5-15) 

where E a is the incident energy in the a-channel and the deuteron wave function 
satisfies 

(e d -h 0a -v a -Lj p )\ip d a ) =0 . (5-16) 
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Note that uj p does not actually contribute in Eq. (j5-16j) . since the isospin of the 
deuteron is zero. From Eq. (|5-15p . we find 

(G p a y X W^) = (E* - h 0a - W p )\ X p a i+) ^ d a ) = . (5-17) 

For three identical particles in the isospin formalism, a transition operator to the 
channel 7, U%, is defined through 

a 

We assume 7 to be the standard coordinate system 7 = 3 and abbreviate the sub- 
script 7. Then, we obtain from Eqs. (|5-10p and (|5T7p the Coulomb-modified AGS 
equation for three identical particles: 

JJP\^+)) = (P V + W P W {+) ) + {Pv + W p )G p U p \ifj p(+) ) , (5-19) 

where G p = (z - H - v - uj p - W p )~ l and P = ^(12)^(23) + f(i3)-P(23) is the 
permutation operator for the rearrangement. In Eq. (|5-19p . we set U p \ip p ^) = 
(Pv + W p )\y p{+) ) and obtain 

= \ippM) + G p {Pv + W p )\V pi+] ) . (5-20) 

Here, is the total wave function for the screened Coulomb problem and is 

related to the total wave function for the full Coulomb problem |<Z^ + )) through PI 

lim W^jC = |^ (+) ) (5-21) 

p— >oo 

with a shift function Q p . The shift function Q p = Q p (k) is defined by 

C p (<Zo) = T W p (R) dR , (5-22) 

where qo is the wave number between the incident proton and the deuteron in the 
center-of-mass (cm) system. The "screening and renormalization procedure"^ con- 
verts Eq. (|5T4p to its full Coulomb correspondence 

{<f>\U c \<t>) = (0|T C |0) + (^H|C/ C |^W) . (5-23) 

Equation (|5-20|) is the distorted-wave version of 

= + g p P( v + u p )\V p(+) ) , (5-24) 



which can be derived similarly from the AGS equation in Eq. (|5-6p by assigning 
U p \4>) = P(v + u} p )\W p ^). In fact, if we note that |^(+)) is three-nucleon anti- 
symmetric, we can easily derive Eq. (|5-19|) from Eq. (|5-24p by using P(v + oj p ) = 



*' Strictly speaking, this relationship is valid only for the sharply cutoff Coulomb potential. 
For general screened Coulomb potentials, an extra finite normalization factor like in Eq. Q4-32p is 
necessary for |!^ p ' + ^). The following relations are all valid by modifying I^M) to \ < P P ^}. 
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Pv + W p + W p . On the other hand, the Faddeev components \ X F P ), satisfying 
|flrP(+)) = (1 + P)\WP), can be derived by setting G U p \<p) = P\V P ) in the AGS 
equation: 

\WP) = |0) + G t p P\V p ) . (5-25) 

In the isospin formalism for the total isospin T = 1/2 state, we use the effective 
T-matrix t I=l = (2/3)i£ p + (l/3)i np for the isospin 1 NN channel.^ 

Instead of using the "screening and renormalization" procedure, we use an exten- 
sion of Vincent and Phatak procedure^ of the two-cluster Coulomb problem, which 
is equivalent to the "screening and renormalization procedure" in the limit of p — > oo. 
The scattering amplitude is obtained by imposing a connection condition on the K- 
matrix0 K p p = {Z" 1 )^ - {(j> a \X p \(j)p) for the pd scattering,® which is derived 
from the two different asymptotic forms of the total wave function in Eqs. (|5-20|) and 
(|5-24p . From here on, the subscripts a, f3, etc. specify the channel quantum num- 
bers. We define a reduced wave function 4>^\r) = (R^ipaWi^) ■ The asymptotic 
form for the wave function Eq (|5-24|) is without a constant normalization factor 

<P p a { p(R) ~ u a (q R) K p n - cv a (q R) 5 aa for R > R out , (5-26) 

where c = q (ir/2)(4M N /3h 2 ) with M N b eing the nucleon mass. For the total wave 
function Eq. (|5-20j) . the asymptotic form is 

$$+\R) ^^J2{ P a(Qo,R)K -cG^ qo ,R)S a , p } 
Qo p 

x — {Fp,i; d \{Pv + W p )\& pi+) ) for R>a, (5-27) 

Qo 

where a is the range of the nuclear force. Here, F p and G p a are the screened Coulomb 
wave functions defined in Eq. (|4-17p . In the inside region R < R[ n , F a and G a are 
equal to F a and G a , respectively. The connection condition for 0^\r) at R = R out 
is written in terms of Wronskians: 

zZKe {w\*$Mro« - W[F$M^ ofay} 

P 

= c{w[G p a ,u a ] Rout K p n - W[G p a ,v a ] Rout c6 a>y ] . (5-28) 

Matrix elements Up , defined by 

'E[KZ t p + ic6 aj f i ]ug„ = 6 a „ (5-29) 
P 

in the limit of p — v oo, are related to (^f ] \U p \^ i+) ) through 

(V$ (-) |fr p |^ (+) ) = e^+^ U p p n . (5-30) 



*' Here, the A"-matrix is defined by the form of Kt(k) = (1/fc) tan 5e(k) for the on-shell matrix 
elements. 
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Here, o~p and cr 7 are the Coulomb phase shifts in the channels j3 and 7, respectively. 

fp'P is obtained from {Jg 



The scattering amplitude is obtained from Uo through 



ha ~ 2 3h 2 to ' [ ' 

In the channel-spin representation, the full scattering amplitude is written as 

t'tJJz 

x ^(£'m'^4| JJ Z ) Y em ,(q f ) Y,(t™S c S cz \JJ z ) Y e * m (q t ) , (5-32) 
m' m 

for a sufficiently large p. 

§6. Numerical performance 

6.1. Comparison with the exact solutions for the Ali-Bodmer aa potential 

Ali-Bodmer aa potential is a simple phenomenological potential which repro- 
duces the results of the phase-shift analysis for the aa scattering up to E cm ~ 15 
MeV. The angular-momentum-dependent version called Ali-Bodmer d (ABd) has 
the explicit form 

V^ d (r) = V l e-^ r2 + V 2 e"^ 2 + —erf (/3r) , (6-1) 

r 

with the parameters 771 = 0.7 2 fm~ 2 , r/2 = 0.475 2 fm -2 , V2 = —130 MeV and 

r 500 MeV for S 
320 MeV for D 
for £ > 4 



Vl 



/o 

8 = - = 0.6014 • • • fm -1 . (6-2) 

H 2 x 1.44 v 1 

In Eq. (|6Tp . erf(x) stands for the error function defined by erf (x) = (2/ y/ir) e - ' 2 dt. 
Since this potential model is exactly solvable by the Runge-Kutta-Gill (RKG) method, 
it is suitable to test the accuracy of the Coulomb approach developed in this paper. 
With the assignment a = 4e 2 , the error function- type Coulomb force 

a 

V D (r) = -erf (Br) (6-3) 
r 

in Eq. (|6Tp is the direct potential of the a-a RGM. When a simple (0s) 4 harmonic- 
oscillator shell-model wave function with the width parameter v is assumed for the 
a-cluster, the parameter j3 is expressed as 



P= 7 -^ = 2^/3, (6-4) 
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where fj, = 4 ■ 4/(4 + 4) = 2 is the reduced mass number of the aa system. On the 
other hand, the rms radius of the a-cluster with A = 4 is given by 



r.-«- 1 |l(l-^-^, (6-5) 
without the proton size effect, so that j3 is related to r a through 

0=277". (6 - 6) 

In ABd, r a = 1.44 fm is assumed, corresponding to v = 0.271 fm -2 . 

In the momentum representation, we use the sharply cut-off Coulomb force at 
the nucleon level. The corresponding direct aa potential is given by 

V&{r) = " | erf (/3r) - i [erf (/3(r + p)) + erf (/3(r - p))] 1 . (6-7) 



If we use this screened Coulomb potential in Eq. (16- ip . we find 

VP a {r) = Vi e-" ir2 + V 2 e~^ r2 + V&[r) . (6-8) 
Here, we separate V£(r) into 

VP{r) = - | [erf {fir) - 1] + 1 - ~ [erf (/3(r + p)) + erf (/3(r - p))\ 



._[l-erf(/3r)] + -a p (r) , (6-9) 

r r 



and set 



a„(r) = l~\ [erf (/3(r + p)) + erf (/3(r - p))] . (6-10) 

Then, the aa potential which should be used in the momentum representation be- 
comes 

V£ a (r) = V(r) + ^a p (r) 

with V(r) = e" r?ir2 + T/ 2 e^ 2 + W(r) . (6-11) 

Here, W(r) = —(a/r) [1 — erf (/3r)] is the short-range attraction originating from the 
Coulomb potential. In fact, the asymptotic expansion of the error function yields 

w^^w-n-VWEt-r^a)" 1 . (6 .i2) 

n=0 V ^ ' 

We find that W(r) is sufficiently small around (f3r) 2 ~ 16; namely, r ~ 4//3 ~ 7 fm. 
(Actually, even around ~ 4 fm, as seen in Fig. [2] below.) 
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Fig. 1. (a): S-wave Ali-Bodmer potential ABd with the screened Coulomb force (solid curve) for 
the aa system. The dashed curve denotes the simple Coulomb potential Vc(r) = 4e 2 /r and the 
dotted curve the nuclear part. The full potential is shown in the bold solid curve. The cut-off 
radius of the sharply cut-off Coulomb force at the nucleon level is assumed to be p = 12 fm. 
(b): The enlarged profiles of (a) for the various Coulomb potentials. The solid curve denotes 
the screened Coulomb direct potential Eq. (|6-7[) with the error function form. 



We illustrate in Fig.[TJa) the S*-wave Ali-Bodmer potential ABd and in (b) the 
enlarged profiles of various types of Coulomb potentials. The cut-off function a p (r) 
in Eq. (|6- 10)) for the cut-off Coulomb radius p = 12 fm and the short-range Coulomb 
potential W(r) in Eq. (|6-12p are shown in Fig.[2j We find that a p (r) satisfies the con- 
ditions 1) - 3) of the screened Coulomb potential. In particular, the much stringent 
condition 3)' in Eq. (|4-5p is also satisfied with the smoothness parameter b ~ 3 fm. 
If we take 6 = 6 fm, the deviation of a p {r) from 1 (or 0) at R\ n = p — b = 6 fm 
(or at -Rout = P + b = 18 fm) is less than 1CP 6 . Note that this kind of a rapid 
transition from 1 to is not achieved in the standard screening functions in the form 
of a p {r) = e~( r / pS) ™ , unless n is taken to be very large like n > 20. In this sense, our 
screened Coulomb potential is a small deviation from the sharply cut-off Coulomb 
potential, which is probably related to the smallness of the limit lim /9 _ >00 ~ 
if it exists. This property must also be related to the small deviation of the shift 
function C, p {k) in Eq. (|4-13p from rylog (2kp), which is the result of the sharply cut-off 
Coulomb potential in Eq. (|2-27[) . We will show in Appendix B that the screening 
function ot p (r) in Eq. (|6-10p satisfies the limit 

C(jfe) -»• V log (2kp) as p^oo, (6-13) 

in contrast to the ce p (r) = e~^ r l case. In the latter case, the right-hand side of Eq. 
(|6-13p contains an extra constant term —(rj/n)j with 7 being the Euler constant. 
(See Eq. ([B^3lU 

First, we have neglected the nuclear potential V\ = and V2 = in Eq. (|6-1|) 
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a p (r)=1-[erf(P(r+p))+erf(P(r-p))]/2 



-1 - 
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Fig. 2. The cut-off function a p (r) in Eq. (|6T0p with the cut-off Coulomb radius p = 12 fm and 
the short-range Coulomb potential W(r)/a = [erf(/3r) — 1] jr in Eq. (|6-12|l for the aa screened 
Coulomb potential (in the unit of fm -1 ). 



and compared the nuclear phase shifts between the present method and the direct 
method using Eq. (|1T|) . In the direct method, the relative wave function ipi(r) in 
Eq. (|l-ip is solved from r = to i? Q ut = 12 + 6 = 18 fm by the RKG method and 
smoothly connected to a linear combination of the pure Coulomb wave functions at 
r = -Rout- Since we are using the error function Coulomb, the nuclear phase shift 
does not become zero, In the S-wave, 5q increases from to 11.088°, when the 
energy increases up to E cm = 15 MeV. Similarly, 5^ = 0.473° and <5f = 0.013° at 
-E-cm = 15 MeV. In the momentum-space approach, we first solve the LS equation 
and calculate 5^ (which is the screened Coulomb phase shift) by assuming p = 12 
fm. The phase shift is then transformed to 5f through the connection condition 
Eq. (|4-24p . Here, we assumed 6 = 6 fm, and F^(k,r), G^(k,r) are calculated from 
R m = 12 — 6 = 6 fm to .Rout = 12 + 6 = 18 fm, also by the RKG method, with the 
pure Coulomb values at R ia = 6 fm as the starting values. The results by these two 
different methods, of course, agree to each other completely within the numerical 
accuracy less than 0.001°. Next, we switch on V\ and V2 and repeated the same 
calculations. The result is shown in Table HI For each incident energy, the first row 
indicates solutions obtained by the RKG method, and the second row those in the 
momentum-space approach. Only different figures from the first row are shown. In 
the left-hand side, the final results of 8¥ are compared. In the right-hand side, the 
phase shifts 5^ directly obtained from the LS equation (before the transformation) 
are also compared. We find that, in the lowest energy E cm = 1 MeV, a difference of 
0.005° exists both in 5f and 8^. This is probably the inaccuracy of solving the LS 
equation in the low energies. For other energies, the difference is less than 0.001°, 
and the agreement of the results obtained by our method with the exact solutions is 
quite satisfactory. 
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Table I. Comparison of the aa nuclear phase shifts {Sf ) of the ABd potential with the direct 
method. For each cm energy E cln , the first row indicates solutions by the RKG method, con- 
nected at i?out = 18 fm by Eq. (|l-ip . The second row stands for the solutions by the present 
momentum-space approach. Only different figures from the first row are shown. In the left-hand 
side, the final results of 5f are compared. In the right-hand side, the phase shifts Sp directly 
obtained from the LS equation (before the transformation) are also compared. The cut-off 
Coulomb radius p is chosen to be p — 12 fm and the smoothness parameter in Eq. (|4-5[) is b = 6 
fm. (h 2 /M N ) = 41.786 MeV • fm 2 and e 2 = 1.44 MeV • fm are used. 
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6.2. aa Lippmann-Schwinger RGM by the Minnesota three-range force 

As a more complex system, we apply the present method to the aa LS-RGM, 
using the Minnesota three-range force. In this calculation, we solve the RGM equa- 
tion in the momentum space. All the Born kernels including the direct term and the 
RGM exchange kernels for the sharply cut-off Coulomb force between two protons 
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Fig. 3. (a): S-, D- and G-wave aa phase shifts predicted by the Ali-Bodmer potential ABd. (b): 
The same as (a), but for aa RGM by Volkov No. 2 (VN2) with m = 0.59 and v = 0.275 fm" 1 
(dashed curves) and by Minnesota three-range (MN3R) potentials with u = 0.94687 and v — 
0.257 fm _1 (solid curves), the latter result is in better agreement with experiment. 



are analytically calculated. For example, the direct Born kernels of the error function 
Coulomb potential in Eq. (|6-3p and the screened Coulomb potential in Eq. (j6-7|) are 
given by 

ME L (Q f ,q l ) = (e^/'H^erf (/3r)|e^-) = 4e 2 ^e-K^) , 

(kp \ 2 \ f h\ 2 
e ~~ 4 ^ } > ( 6 - 14 ) 

where k = qf — q t . Note that Mp L (qj, involves the Coulomb singularity at \q^\ = 

\qi\-, while L (qf,qi) does not have such a singularity. A numerical challenge is 
the angular momentum projection of this kernel. We have used a standard Gauss- 
Legendre integration quadrature, taking many discretization points. We can check 
the accuracy of this numerical integration by examining the redundancy condition 
of the Pauli forbidden states for the S- and L>-waves. Various cut-off Coulomb 
parameters are chosen from p = 8 fm to 16 fm, with b = 6 f m fixed, The modified 
Coulomb wave functions are therefore solved from R[ n = p — 6 fm to R ut = p + 6 fm. 
In Table 2, we list the variation of the nuclear phase shifts, depending on the choice 
of p. We find that the results are quite stable in this appropriate range of p. We 
show in Fig.[3^a) the aa phase shifts predicted by Ali-Bodmer d potential and in 
Fig. 0(b) the results by the LS-RGM using the Minnesota three-range force and the 
Volkov No. 2 two-range force. 

6.3. pd elastic scattering 

As in the case of the aa scattering discussed in the preceding subsections, the 
screening function a p (R) for the pd elastic scattering should be derived in a consistent 
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Table II. Cut-off radius (p) dependence of the nuclear phase shifts 5e for the act Lippmann- 
Schwinger RGM. The Minnesota three-range force with u — 0.94687 and v = 0.257 fm -2 
are used. 
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4.000 


0.398 


0.399 


0.399 


0.399 


0.399 




5.000 


1.107 


1.108 


1.108 


1.108 


1.108 




6.000 


2.522 


2.523 


2.523 


2.523 


2.523 




8.000 


9.514 


9.517 


9.517 


9.517 


9.517 




10.000 


29.924 


29.928 


29.928 


29.928 


29.928 




12.000 


75.818 


75.820 


75.820 


75.820 


75.820 




15.000 


120.744 


120.745 


120.745 


120.745 


120.745 



way with the screened Coulomb potential between two protons in Eq. (|5-ip . In 
our application of the quark-model baryon-baryon interaction fss2 to the pd elastic 
scattering in Ref. 29), the sharply cut-off Coulomb force is introduced at the quark 
level in the form of {l/r qq )0(p — r qq ), where r qq is the relative distance between two 
quarks. The proton-proton (pp) potential u p (r) is obtained by folding it with the 
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(3q)-(3q) internal wave function, resulting in 



w'(r;l,2) 



e 
r 



erf 



V3 r 
~2~b 



erf ( — ] + erf 



y/3: 



1 + r z (l) 1 + r z (2) 



(6-15) 



where r is the distance between the two protons, r = |-r 12 1 = | £Ci — a^l and b is the 
harmonic-oscillator range parameter of the (3(7)-clusters. Note that this screened 
Coulomb potential for the two protons is not equal to Eq. (|5Tp with a mere change 
of a p (r) to Eq. (|6TU|) (with a trivial modification f3 — > (\/3/26)), but also contains 
the contributions from the short-range Coulomb potential in Eq. (|6T2p . We calculate 
the pd screened Coulomb potential by further folding the pp potential in Eq. (J6T5P 
with the deuteron wave function (r; 1, 2\ipd)'- 



V p f(R) = (r\^(\R + r/2\;2,3)\r) + (r\^(\R-r/2\;3,m d ) , (6-16) 

where R = X3 — (x\ + x<i)j2 is the relative coordinate between the center-of-mass 
of the deuteron and the proton. This calculation is made in Appendix C. We assign 
the long-range part of Vf^'(R) in Eq. (|C-2jl to W P (R) in Eq. (|5-2p . and parametrize 
it as W P (R) = (e 2 / R)a p (R). The screening function a p (R) is numerically calculated 
by using Eqs. (|CT5p - (|CT7p and the momentum-space deuteron wave function ex- 
panded in the dipole form factors!^* Here, we only show in Fig. 2] the profiles of the 
screening function a p (R) and the short-range Coulomb potential (the polarization 
potential) W(R) for the simplest deuteron channel with J 71 = l/2 + . We find that 
the cut-off behavior around R ~ p is fairly sharp even in p ~ 8 fm. The short-range 
Coulomb potential W(R) is ^-independent as shown in Eq. (|(J-5|) . The coupling 
potential W P (R) between different channel-spin states, (£S C ) 7^ (£' 'S' c ), is very small. 
We therefore neglect this and solve the screened Coulomb problem only by using 
the diagonal part of (£S C ), in order to generate the regular and irregular screened 
Coulomb wave functions for the connection condition. 

Some typical eigenphase shifts of the E p = 65 MeV pd scattering with the 
Coulomb cut-off radius p = 8, 16 and 20 fm are listed in Table IIIII for the S and 
P waves. Here, we have assumed the maximum total angular-momentum of the 
two-nucleon subsystem, / max = 4. The real parts of the eigenphase shifts are only 
given for simplicity. We find that the inclusion of the cut-off Coulomb force gives 
an apparent repulsive effect, namely, the S- wave and P-wave eigenphase shifts are 
—0.9° ~ —2.5° (—2.5° ~ —3.1°) more repulsive than in the no Coulomb case, if 
p = 8 fm (p = 16 fm) is assumed. The transformation by the connection condition 
in Eq. (|5-28p gives an attractive effect to make the resultant eigenphase shifts rather 
close to the no Coulomb case. As long as the low partial waves such as the S and 
P waves are concerned, the final results of the nuclear eigenphase shifts are rather 
stable within the fluctuation of less than 0.8°. We have calculated pd differential 
cross sections and other polarization observables, using various p values. The results 
by p = 8 fm is quite reasonable, but if we take larger values like p = 16 and 
20 fm, we have found that unpleasant oscillations develop in all the observables. 
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4. (a) The screening function a p (R) (solid curve) and the short-range potential W(i?)/e 2 
(dashed curve) (in the unit of fm -1 ), given in Eq. (|C-7[I for the pd scattering. A simple 5- wave 
deuteron wave function u(r) = \Z2j e ,r is used. The Coulomb cut-off radius is p = 8 fm. (b) 
The realistic a p (R) for the pd scattering in the simplest deuteron channel with J" = l/2 + . 
In the (£S C ) = (0 1/2) (solid curve) and (2 3/2) (dashed curve) diagonal channels, curves are 
almost identical. The off-diagonal a p (R) with (iS c )-(i'S' a )=(0 l/2)-(2 3/2) is very small. Here, 
£ is the relative angular-momentum between p and d and S c is the channel spin. 



The origin of the oscillations is traced back to the high partial waves, in which the 
restriction of / max = 4 is too severe. Since we are using the channel-spin formalism, 



Table III. Real parts of the nuclear eigenphase shifts for the Nd elastic scattering at Em = 65 MeV. 
The nd phase shifts with no Coulomb force and the pd phase shifts including the cut-off Coulomb 
force with p = 8, 16 and 20 fm are listed. For p = 8 fm (before), the eigenphase shifts before 
transformation in Eq. (|5-28p are also shown. The maximum total angular-momentum of the two- 
nucleon subsystem is I m&x = 4, and the momentum discretization points n — n\-ni-nz = 6-6-5 
are used in the definition shown in Ref.l28[). 



2S+l £j 


no Coulomb 


with Coulomb 


p = 8 fm 
(before) 


p = 8 fm 


p = 16 fm 


p = 20 fm 


2 Si/2 


26.84 


24.38 


28.70 


28.99 


29.01 


*D 1/2 


-7.25 


-9.41 


-7.20 


-6.92 


-6.76 


2 p 

-fl/2 


-0.44 


-2.51 


-0.04 


0.35 


0.45 


" 1/2 


24.28 


21.83 


24.76 


24.98 


24.99 


4 'S'3/2 


32.11 


31.23 


33.79 


34.25 


34.63 


2 D 3/2 


8.74 


6.87 


9.15 


9.35 


9.51 


4 D 3/2 


-5.49 


-7.84 


-5.32 


-5.01 


-5.11 


4 P3/2 


24.98 


22.53 


25.32 


25.84 


25.72 


2 ^3/2 


6.73 


4.60 


7.16 


7.38 


7.47 


4 F 3/2 


-1.05 


-2.81 


-0.86 


-0.63 


-0.55 



36 



MS) 



Y. Fujiwara and K. Fukukawa 
A y (6) 



0.06 



0.04 



0.02 



E p =3 MeV 
n=6-6-5 

I -3 
'max 



dot-dot-dash no Coulomb 

dot 8 fm 
dash 10 fm 

o?° 

-9 solid12fm 

a< — . B 
oh 




0.06 



0.04 



0.02 



180 







E =3 MeV 


dot-dot 


-dash no Coulomb 


n=6-6-5 
I =4 

'max 

O' 

"// 

0, / 

o°// 
°o,'/ 

CO g 


dot 8fm 
dash 1 fm 
oo solid 12 fm 

\\° 

V 

V 


o 






\ vO 








/'S i 











60 



120 



180 



e C m(deg) e cm (deg) 

Fig. 5. Cut-ofT radius dependence of the proton analyzing power for the pd elastic scattering at 
E p = 3 MeV. The results of no Coulomb case, p = 8, 10, and 12 fm are shown by the dot- 
dot-dashed, dotted, dashed, and solid curves, respectively. The left panel shows the results of 
/max = 3 and the right panel of / max = 4. The pd experimental data from Ref.l32p are also 
shown by circles. 



the total angular momentum J 71 " of the three nucleon system is achieved by the 
angular-momentum coupling (£S C )J, where the channel spin S c is constructed from 
(I^)S C . For a large J 77 , a large contribution of the Coulomb force from the large 
relative orbital angular momentum of the two-proton subsystem is not fully taken 
into account, since the magnitude of S c is restricted by I max = 4. To demonstrate 
this situation, we show in Fig. [5] the /^-dependence of the nucleon analyzing power 
for the 3 MeV pd scattering, calculated with J max = 3 and I max = 4. In the forward 
angular region with 9 cm < 90° , we find that unpleasant bump structure develops as p 
increases from 8 fm to 12 fm, when / max = 3 is used. However, such enhancement is 
strongly suppressed when J max = 4 is used. This demonstrates very clearly that two- 
nucleon partial waves should be included up to sufficiently higher values to obtain 
the well converged results, if the screened Coulomb force is incorporated into the 
standard AGS equations. 

Since the calculation with I max = 6 and more is not presently possible because 
of the computer resources, here we propose to cut the Coulomb force for higher J w 
values and use a simple "Coulomb externally corrected approximation" , in which the 
nd eigenphase shifts are directly used for the nuclear phase shiftsPD Figure \U\ shows 
the pd differential cross sections and some polarization observables at E p = 65 MeV, 
calculated by neglecting the Coulomb force for J n > ll/2 + . We find that the results 
with p = 8, 16 and 20 fm are very similar, although some difference is seen in T22{6)- 
The results with p = 8 fm are almost the same as the full calculation including the 
Coulomb force to all the partial waves. 
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Fig. 6. pd differential cross sections (do/dfi), analyzing power (A y (6)) of the proton, and vector- 
type (iTn(8)) and tensor-type (T2 m (f?)) analyzing powers of the deuteron at E v — 65 MeV. 
The results with no Coulomb case, p = 8, 16, and 20fm are shown by the dashed, dotted, 
solid, and bold-solid curves, respectively. These curves almost overlap with each other, except 
for the forward nuclear-Coulomb interference region. The screened Coulomb force is neglected 
for higher partial waves with J w > ll/2 + . The experimental data are taken from Ref. l33p for 
dajdQ and A v (6), and from Ref.|34| for iTu{9) and T 2m (6»). 
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§7. Summary and outlook 

In the present work we have proposed a practical method to deal with the 
Coulomb problem in the momentum space. Although the standard procedure to deal 
with the Coulomb force in two-body systems is formulated in the configuration space, 
the extension of such an approach to three-body systems is not trivial.® Here, we 
have reformulated the momentum-space approach of the two-cluster systems based 
on the essential idea of the "screening and renormalization procedure" , which is re- 
cently used in the standard formulation of the AGS equations for the pd scattering in 
the momentum representation.®'®'® In this approach, the screened Coulomb force 
with a cut-off parameter p is introduced to the basic equations as if it is a part of the 
short-range nuclear force. The two-potential formula for the short-range potentials 
is used to generate the scattering amplitude. The pure Coulomb results are repro- 
duced by taking p — > oo limit, based on the Taylor's formula^® for the phase 
renormalization of the asymptotic wave functions of the screened Coulomb poten- 
tial. The central issue in this approach is if one can reproduce the exact Coulomb 
results by taking a finite p. Since the quasi-singular nature of the screened Coulomb 
force becomes stronger for larger p, it is essential that one can reproduce the almost 
exact results with a reasonable choice of p. 

To achieve this, we propose to extend the Vincent and Phatak approach, 6 ' which 
is originally formulated for the sharply cut-off Coulomb problems. When a sharply 
cut-off Coulomb force with the cut-off radius p is introduced at the level of constituent 
particles, two-cluster direct potential of the Coulomb force becomes in general a local 
screened Coulomb potential implemented with the short-range Coulomb force. The 
screening function a p (r) is determined by the properties of the cluster wave functions, 
and involves a smoothness parameter b related to the size of clusters. In practice, 
b satisfies b -C p, which is an additional condition to the Taylor's properties^® of 
screening functions. We find that this condition is necessary to make the present 
treatment work well. We pay attention to the existence of two different types of 
asymptotic waves contained in the screened Coulomb wave functions. The first one 
is the approximate Coulomb wave for the relative distance of two clusters, r, smaller 
than i?j n = p — b, and the other is the free (no-Coulomb) wave in the longer range 
region, r > -Rout = p + b. The asymptotic Hamiltonian composed of the screened 
Coulomb force allows us to calculate constant Wronskians of this Hamiltonian in 
either region. Using this property, we can extend the standard procedure of matching 
conditions for asymptotic waves to the screened Coulomb potential. 

We should note that the renormalization property of the screened Coulomb wave 
functions is more involved than in the sharply cut-off Coulomb case. In particular, 
the irregular function of the screened Coulomb potential in general contains an ad- 
mixture of the regular solution even in the p — > oo limit. As the result, the limit of 
the Green function for the screened Coulomb potential is not reduced to the Coulomb 
Green function. This requires an extra renormalization of the regular wave function 
for the problem of the short-range nuclear potential plus the screened Coulomb po- 
tential. This renormalization factor, however, does not affect the final expression of 
the connection condition, since it is given by the ratio of Wronskians. 
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We have applied this method first to an exactly solvable model of the aa scat- 
tering with the Ali-Bodmer potential and confirmed that essentially exact phase 
shifts are reproduced, using a finite p. The stability of nuclear phase shifts with re- 
spect to the change of p in some appropriate range is demonstrated by using the aa 
Lippmann-Schwinger RGM with the Minnesota three-range force. In the application 
to the pd elastic scattering, some dependence on the choice of p remains although 
the essential features of the nuclear and Coulomb interference in forward angles are 
reproduced not only for the differential cross sections but also for the deuteron tensor 
analyzing powers. 

We have to admit that the completely satisfactory Coulomb treatment of the 
three-body system is still beyond the way. First, the stability of p in the case of 
the above pd elastic scattering is not completely realized. We have examined all the 
observables for the pd elastic scattering in the energy range E p < 65 MeV, and found 
that the present choice p ~ 8 - 9 fm is a reasonable choice to reproduce almost all 
the experimental data.® The forward behavior of the vector analyzing power A y (9) 
for the proton and iT\\{9) for the deuteron is not consistently achieved in the low- 
energy region, using a unique p. Choosing much larger p around p ~ 16 - 20 fm is 
almost prohibited since the solution of the AGS equation becomes very singular and 
the partial waves included in the actual calculations are restricted by the hardware. 
Another problem is the treatment of the Coulomb force in the breakup processes.^ 1 
The phase renormalization for the observed two protons at the final stage is not 
trivial because of the exchange breakup amplitude. We probably need to solve the 
Coulomb-modified AGS equations in spite of the very singular nature of the screened 
Coulomb wave functions in the momentum representation. Finally, we mention that 
the Coulomb treatment of three charged particles like the three a system is a big 
challenge, since the the asymptotic behavior of the three charged particles is not 
a prior known. 
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Appendix A 

Definition of the Coulomb wave functions 

The usual regular solution ipe(k,r) and the regular solution corresponding to 
the Jost solution, (pg(k,r), for the Coulomb problem are defined by the confluent 
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hypergeometrical functions through 



,£+1 

-e tkr F{£ + l+i V ,2£ + 2, -2ikr) 



(2£ + l)U 

= (2l + l)\\ e ~ lkrF{f ' + 1 -^ 2£ + 2 > 2ikr ) = real 
~ ynp[\ F i{k)\ sin (kr - rj\og2kr - ^£ + a^j (r oo) . (A-l) 
Here, n is the Sommerfeld parameter and i^(fc) is the Coulomb Jost function 

F e (k) = 11 , (A-2) 

w r(£+l + ir)) ' v ' 

which can be obtained by comparing the behavior at the origin between ^i(k,r) 
and <pi(k,r). The Jost solution of the Coulomb problem is defined by the irregular 
solution with the asymptotic behavior fg(k, r) ~ e i{kr-r)\og2kr-iri/2) £ Qr r _^ More 
explicitly, it is given by 

fi(k, r) = {-if{2kr)-^e lkr G{£ + 1 + irj, -£ + ir], 2ikr) 

= i(-Y +1 e nv/2 {2kr) i+1 e ikr ^(£ + l + ir),2£ + 2, -2ikr) 

„ el (kr- V log2kr~ne/2) ^ _^ ^ _ ^.3) 

Here, G(a, f3, z) and #(a, 7, z) are irregular solutions of the confluent hyper geometric 
functions defined in Refs.122]) and [23]), respectively, and they are related to each other 
by 

#(0,7,2) = z~ a G(a,a -7 + l,-z) . (A-4) 
The symmetries of the Jost solution and the Jost function are given by 

f* £ (k,r) = (-)' e ™>f t (-k,r) , F/(fc) = e™>F t (-k) , (A-5) 
with the Coulomb factor e 7 ™ 7 . They satisfy the usual definition of the Jost function 



and the relationship 



Mk,r) = ^^{FU^mr) - Fe^fUKr)} , (A-7) 

for real k. 

The usual Coulomb wave functions are defined as the real functions satisfying 
the asymptotic behavior 

Fg(k, r) ~ sin (kr — rjlog 2kr — —£ + a^j , 

Ge(k, r) ~ cos [kr - rjiog 2kr + 01) > (A-8) 
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for r — > oo. These Coulomb wave functions are related to each other through 



ip £ (k,r) = ~e ia *F e (k,r) 



fe(k,r) 



k i+i 



\F e (k)\ F e (k,r) =real , 
G e (k,r) + iF e (k,r)] , 



ft(k,r) = e i ^[G e (k,r)-iF e (k,r)] 



(A- 



The relationship with the usual "incident plane wave + outgoing (or incoming) 
spherical wave" is ipp~\k,r) = ipt(k,r) and tpi \k,r) = ip^(k,r). This implies that 

(k, r) = tpe{k, r) ~ 'i.^ 1 s ™ ~ V^°E %kr — —£ + a^j , 



7T 



pC i(kr--q\og2kr~(-K /2)i) 



~ — sin I kr — r]log 2kr 1 ) + fp e 

k \ 2 / 

ip£~\k, r) = tpe (k, r) ~ \e~ i(Tl sin (kr - rjlog 2kr -~£ + ae 



pC* -i(fcr-»jlog 2kr-{-n/2)l) 



~ — sin I fcr — rylog 2/cr 1 ) + fg e 

Here, /f = (l/2ife)(e 2i ^ - 1) is the Coulomb partial -wave amplitude, and 



(A40) 



D 2ia e 



r(l + 1 + i V ) 
r(£ + l- irj) 



\m)\ 



,2wri _ j 
2tXT] 



n 



^ n 2 + r/ 2 



(A-ll) 



Appendix B 

57u/i function of various screening functions 



In this appendix, we calculate the shift function 



C p (k) 



1 

2k- 



2kr] 



-a. 



i_ r 

2k 



(r) dr = r/ —a p (r) dr , 



2fe 



(B-l) 



appearing in Eq. (jl -5j) for various screening functions a p (r) and evaluate the no- 
screening limit p — > oo. When the screening is a p (r) = e - ( r / p )", we can write an 
analytic expression 

coo 



C(k) = r, 
which leads to 



2k 



e"(p) dr = r) log (2fcp) - ? 7 + - ^ f - 1 



n n 



~i r r\ \2kp 



(B-2) 



C p (&) — >■ ?7 log (2fc/o) 7 as p — > oo 

n 



(B-3) 



Here, 7 is the Euler constant. On the other hand, the screening functions with more 
sharp transitions like 3)' in Eq. (|4-5p seem to have no constant term like Eq. (|6-13|) 
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in the limit of p — > oo. We will show this for the error function screening in Eq. 
(|6-10p . The proof for the exponential screening function in Eq. (|C-7|) is also carried 
out similarly. 

In order to prove Eq. (|6- 13j) . we separate the r integral in Eq. (|B-ip into three 
pieces as 

fP 1 fP 1 f°° 1 

( p (k) = rj - dr — 77 / -(1 — a„(r)) dr + rj / -a p (r) dr 

= 77log(2^)-/i(p)+/ 2 (p) • (B-4) 
First, the positive integral hip) is estimated by 

/ 2 (p) < - / dr , (B-5) 

so that we only need to evaluate the integral over a p (r). For the error function 
screening, the expression 



a p {r) = -^={ I e' tz dt + / e~*" dt }> (B-6) 



TT 



yields 



a p (r)dr = - 7 =-(l + e-( 2 ^ 2 )-|(l-erf(2/3 /3 )) . (B-7) 



We therefore find 

/2(p)< ^2^( 1 + e " (2/3P)2 )"I (1 " erf(2/3/3)) ^° aS P ^°°- (B ' 



ttp (r) = 1 - / e"* 2 (it (B- 

/3(p-r) 



In order to evaluate h(p), we use 

1 /-j8(p+r) 

A 

derived from Eq. ()B-6j) . and express it as 

Ii(p) = 4= / - / e" 1 dt dr. (B-10) 

Here, we change the integral variable from r to x by r = /jx and obtain 

n f 1 1 / /•Q!(l+a:) \ 

J x (p) = -L / - / e ~* dt dx , (B-ll) 

with a = fip and e = We consider the upper bound ^Ji(a) > Ji(p) with 
~ f 1 1 / \ 

ii(a) = / - / e _t dt dx . (B-12) 

JO x lia(l-i) / 
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We separate the integral interval [0, 1] into [0, 1 — 5] and [1 — 5, 1] with a small positive 
5 > 0. Then, we find 

1 / ra(l+x) \ 
I^a) = - / e~* dt ] dx 

JO x \Ja(l-x) J 

+ / - / e~* dt) dx . (B43) 

Here, the first term is bounded by 2ae~* a5 * )2 (1 — 5). As to the second term, we 
change the integral variable from x to y by x = 1 - y and find 

2-nd term = / / e~* dt] dy 



1-1/ 



at/ 



Thus, we obtain 



< Ix(a) < 2ae-^ 2 (l - 5) + 

l — o 2 

First, we take the limit a — > oo in Eq. (|B-15p and obtain 



(B45) 



< lim h(a) < ^— x ^- . (B-16) 

a^oo 1 — 2 



Since we can take 5 > arbitrary small, we eventually find 



h(p) < -^=Ii(a) — ► as oo . (B-17) 

\/7T 



Appendix C 

Screening function a p (R) for the dp scattering 

In this appendix, we derive the screening function a p (R) for the pd scattering, 
starting from the screened Coulomb function w p (r;l,2) in Eq. (|6T5p for the pp 
system of the quark-model baryon-baryon interaction. We first note that the p — > oo 
limit, uj = linip-j.oo lo p , is an error function Coulomb, which satisfies 

<^|(Pw)|^) ~ ^ for R -> oo , (C-l) 

where (1,2|^) is the deuteron wave function and P is the rearrangement permuta- 
tion operator P = P(i2)-P(23) + -P(i3)-P(23)- We follow the procedure similar to the aa 
case in Eq. (j6-9|) and separate the folding pd potential in Eq. f)616f) for the screened 
Coulomb force into the long-range and short-range parts: 
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e 2 ,„ _ , A / e 2 



- - (MPw) " (^ P )I^)J + ( + (^l(^)l^) 
e 2 

= — a p ( J R)+W( J R) = P^ p ( J R) + W( J R) . (C-2) 

Here, the screening function a p (R) and the short-range Coulomb potential W(R) is 
given by 

a?>(R) = l-*(il> d \(Pu)-{Pu p )m , 

W{R) = m{Pu)- e ^ d ) ■ (C-3) 



On the other hand, the exchange term Eq. (|5-2p in the three-body model space yields 
the matrix element 

v P P d( R ) = &d\W p + W\il> d ) = (VdlWlVd) + WP{R) . (C-4) 

We therefore find that the deuteron matrix element of the polarization potential is 
p- independent: 

WdWWd) = w(R) • (c-5) 

We first assume the sharply cut-off Coulomb force Eq. (|5-ip with a p (r) = 9{p—r) 
for the two protons, and examine the screening property discussed in §4 by using 
available analytic expressions. This is possible, if we further neglect the .D-state 
component of the deuteron wave function and assume that the spatial part of the S- 
wave component is given by a simple exponential function u(r) = y / 2 7 e~ yr . In this 
case, the folding potential is expressed in terms of the integral exponential function 
defined by 

f°° e ~ f 
Ei(-x) = - / dt 

Jx t 

x 2 ( — x) r 

= log x + 7 - x + — — h -. < (for x > 0) 

2-2! r ■ r! 

n=l 

We find 



~e x ^2(—) n — (asymptotic expansion) . (C-6) 



e 2 e 2 



W(R) = V p c d (R) - C - = -\e~^ R - 4 7 e 2 Ei(-4 7j R) 



-2 



-e" 47i? — ^ I 1 — 7—=: + 



2! 3! 



4 7 i? 2 V 4 7j R (4 7j R) 2 
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poo 

a p (R) = 1 - 2 7 / dr [ Ei(-4 7 (r + R)) - Ei(-4 7 |r - R\) ] 

J O 
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-27 



rcxj rcxj 

/ dr Ei(-4 7 r) + / dr Ei( 
./fl-p JR+p 



-4 7 r) 



for i? > p . (C-7) 



Here, V^(i?) = linip^oo V£ d (R) = (ipd\(Pui)\i/j,i) . In order to derive the last expres- 
sion of Eq. (|C-7j) , we use the relationship 

pOO 

(-4 7 ) / dr Ei(-4 7 r) = 1 , (C-8) 
J o 

which is obtained by exchanging the integration order. From here, we can obtain 

poo poo 

(-27)/ dr Ei(-47|r-J?|) = 1 + 2 7 / dr Ei(-4 7 r) for i? > p . (C-9) 



The asymptotic form of W(R) in Eq. ()C-7j) is due to 



-4 7 ) Ei(-4 7 r) ~ -e~ 47r ( 1 - — + 



2! 



47r (47r) 2 

If we further use this in the last expression of Eq. (|C-7p , we find 



<*p(R) ~ " g t Ei(-4 7 (i? - p)) + Ei(-4 7 (i? + p)) 



as r — > 00 . (C-10) 



-Ei(-4 7j R) ^ as i?^oo 



(C-ll) 



With R fixed, we can show linip^oo a p (R) = 1 as follows. First, Eq. ()C-9p and some 
calculations yield 



27 r°° 

dr [ Ei(-4 7 (r + .R)) - Ei(-4 7 |r - ] 

R Jo 

= ^ + ^f ° dr Ei(-4 7 r) = ^ (1 - e" 4 ^) - 4 7 Ei(-4 7 i?) . (C-12) 
Thus, Eq. ([UTOD gives 

poo 

2 7 / dr [Ei(-4 7 (r + i?))-Ei(-47|r-i?|) ] 



l-e 



1 



1 



2! 



4 7 i? V 4 7-R 



as R — > 00 



(CMS) 



Here, because of r + R > \r — R\, the integrand of Eq. (|C-13P is always positive. 
Furthermore, the integral from to p in Eq. (|C-13p is the monotonically increasing 
function of p and the limit p — > 00 exists. We therefore find 

poo 

lim 2 7 / dr [ Ei(-4 7 (r + i?)) - Ei(-4 7 |r - ] = . (C-14) 
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After all, we find that the screening function a p {R) satisfies the condition 1) - 3) at 
least in this simplest case. If we calculate the shift function ( p (k) using a p {R) in 
Eq. (|C-7p . we find the same result Eq. (|6-13|) ; namely, there is no constant term as 
in the sharply cut-off Coulomb case. 

The calculation of a p {R) using the screened Coulomb potential in Eq. (|6- 15|) and 
the realistic deuteron wave function by the quark-model baryon-baryon interaction 
is rather involved. We here show only the final result for the numerical calculations. 
The screening function a p {R) in the channel-spin formalism is given by 



A,A'=0,2k=0,2,4 



(C-15) 



where the kinematical factor 9 S') * s gi ven by 

9(£SaW's' c ) = {-) Sc+s ' c+1 ^sJ' c x^'e'(^ox'o\KO)(me'o\HiO) 
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J £' S' c 
X' S L 



(C-16) 



Th spatial function (R) is given by 



fZ x ,(R) = dr u x (r)uy(r)v K (R,r/2) , 
J o 

v K (R,r/2) = ^J dx v(^/R 2 + r 2 /A-rRx} P K (x) , 



v(r) 



R 
2r 



erf 



V3r + 



+ erf 



y/3; 



(C-17) 



where u\(r) is the S'-wave (A = 0)) and .D-wave (A 
usually denoted by u(r) and w(r), respectively. 



2) deuteron wave functions 
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